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ABSTRACT 

We have measured gluon splitting into bottom quarks, g — ► bb, in hadronic Z° decays collected 
by SLD between 1996 and 1998. Lhe analysis was performed by looking for secondary 
bottom production in 4-jet events of any primary flavor. 4-jet events were identified, and 
in each event a topological vertex-mass technique was applied to the two jets closest in 
angle in order to identify them as b or b jets. The upgraded CCD-based vertex detector 
gives very high i?-tagging efficiency, especially for B hadrons with the low energies typical 
of this process. We measured the rate of g — > bb production per hadronic event, g b ^, to be 
(2.44 ± 0.59(stat.) ± 0.34(syst.)) x 1(T 3 . 
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1. Introduction 



The vertex representing a gluon splitting into a heavy-quark pair, g — > QQ {Q = b or 
c), is a fundamental elementary component of Quantum Chromodynamics (QCD), but the 
contribution of this vertex to physical processes is poorly known, both theoretically and 
experimentally. In high-energy e + e~ annihilation the leading-order process containing this 
vertex is e + e~ — ► qqg — > qq QQ (q = u, d, s, c or b). Information on g — > QQ can hence be 
obtained by studying e + e _ — > hadrons events comprising four quark jets, with two of the 
jets identified as Q or Q. Background events of the kind e + e~ — > QQ g — > QQ are in 
principle indistinguishable final states, although their kinematics are typically quite different: 
the Q and Q jets in this process tend to be back-to-back, and have energy comparable with 
the beam energy. By contrast, the Q and Q jets from g — ► QQ tend to be collinear, and 
have low energy. In order to enhance the contribution from the latter process it is necessary 
to identify four-jet events and to tag typically the two jets closest in angle, and/or of lowest 
energy, as Q or Q. 

We define the rate gnn as the fraction of e + e _ — > hadrons events in which a gluon splits 
into QQ, e + e~ — > qqg — > qq QQ. Since the quark mass provides a natural cutoff gqq is an 
infrared finite quantity, which can be computed in the framework of perturbative QCD. At 
the Z° resonance energy the production of secondary cc or bb via gluon splitting is strongly 
suppressed by the required large gluon virtual mass. From a leading-order + next-to- leading- 
logarithm approximation calculation one expects |2j g cS to be at the per cent level, and g b i 
to be at the 0.1 per cent level. The calculations, however, depend on a s and on the quark 
mass, which results in substantial theoretical uncertainties. 

Here we consider the process g — ► bb. The measurement of g b i is difficult experimentally 
since the rate is intrinsically low and the backgrounds from Z° — ► bb events are two orders 
of magnitude larger. In addition, the B hadrons from g — > bb have relatively low energy 
and short flight distance and are difficult to identify using standard tagging techniques. So 
far, measurements of g hb have been reported by ALEPH, DELPHI and OPAL ||; the most 
precise measurement, from ALEPH, is g bb = (2.77 ± 0.71) x 10 -3 . 

Such limited knowledge of g bb results in the main source of uncertainty in the measurement 
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of the partial decay width R b = T(Z° — > bb)/T(Z° — > qq) Q, which is potentially sensitive, 
via loop effects, to new physics processes that couple to the fo-quark. Hence, more precise 
measurements of g bb would help to improve the precision of tests of the electroweak theory 
in the heavy-quark sector. In addition, knowledge of the g — > bb process is vital for mea- 
surements at hadron colliders. For example, about 50% of the B hadrons produced in QCD 
processes at the Tevatron are due to g — > bb, and a larger fraction is expected to contribute 
at the LHC. These events form a large background to possible rare new processes involving 
decays to heavy quarks, such as H° — > bb, and improved understanding of g — > bb will help 
to constrain background heavy-flavor production at hadron colliders. 

We present a measurement of g bb based on the sample of roughly 400,000 hadronic Z° 
decays produced in e + e~ annihilations at the Stanford Linear Collider (SLC) between 1996 
and 1998 and collected in the SLC Large Detector (SLD). In this period, Z° decays were 
collected with an upgraded vertex detector with wide acceptance and excellent impact pa- 
rameter resolution, thus improving considerably our tagging capability for the low-energy B 
hadrons characteristic of the g — > bb process. 

2. The SLD 

A description of the SLD is given elsewhere . Only the details most relevant to this anal- 
ysis are mentioned here. The trigger and selection criteria for Z° — > hadrons events are 
described elsewhere ||. This analysis used charged tracks measured in the Central Drift 
Chamber (CDC) @ and in the upgraded CCD Vertex Detector (VXD) ||, with a momen- 
tum resolution of cr p± /p± = 0.01 © 0.0026pj_, where p± is the track transverse momentum 
with respect to the beamline, in GeV/c. For high-momentum tracks the measured impact- 
parameter resolution approaches 7.7/j,m (9.6//m) in the plane transverse to (containing) the 
beamline, while multiple scattering contributions are 29//m/(psin 3 / 2 9) in both projections. 
In e + e~ — > hadrons events the centroid of the SLC interaction point (IP) was reconstructed 
with a precision of approximately 5/im (10/xm). 

Only well-reconstructed tracks |J] were used for i?-hadron tagging, and tracks from iden- 
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tified 7 conversions and K° or A decays were removed from consideration. Each track was 
required to have: a polar angle satisfying |cos#| < 0.87, a transverse impact parameter < 
0.30cm with an error < 250/im, impact parameters, measured in the CDC only, of < 1.0cm 
(transverse) and < 1.5cm (plane containing the beamline), at least 23 hits in the CDC with 
the first hit < 50.0cm from the IP, a x 2 /d.o.f. < 8.0 for the CDC-only and CDC+VXD fits, 
and p L > 0.25GeV/c. 

For the purpose of estimating the efficiency and purity of the g —* bb selection proce- 
dure, we made use of a detailed Monte-Carlo simulation of the detector. The JETSET 7.4 



[ 10 1 event generator was used, with parameter values tuned to hadronic e + e~ annihilation 
data [11], combined with a simulation of B hadron decays tuned to T(4S*) data |]l2[ and a 
simulation of the SLD based on GEANT 3.21 IfTHfl . Inclusive distributions of single-particle 
and event-topology observables in hadronic events were found to be well described by the 
simulations ||14| . Uncertainties in the simulation were taken into account in the systematic 
errors (Section]^). 

Monte-Carlo events were reweighted to take into account the current measurements of 
gluon splitting into heavy-quark pairs H, |15|| . JETSET with the SLD parameters predicts 
g b l = 0.14% and g c5 = 1.36%. We reweighted these events in the simulated sample to obtain 
g bi = 0.247% and g c5 = 3.07% @. Samples of about 1900k Monte-Carlo Z° -» qq events, 
1900k Z° — > bb events, 1090k Z° — > cc events and 60k Z° — > qqg, g — > bb events were used 
to evaluate the selection efficiencies (Section ^). 



3. Flavor Tagging 



We used topologically-reconstructed secondary vertices [17] for heavy-quark tagging. To 
reconstruct the secondary vertices, the space points where track density functions overlap 
were found in three dimensions. Only the vertices that are significantly displaced from the 
IP were considered to be possible B- or D-hadron decay vertices. The mass of the secondary 
vertex was calculated using the tracks that were associated with the vertex. We corrected the 
reconstructed mass to account for neutral decay products and tracks missed from the vertex. 
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By using kinematic information from the vertex flight path and the momentum sum of the 
tracks associated with the secondary vertex, we calculated the Py-corrected mass, Mp T , by 
adding a component of missing momentum to the invariant mass, as follows: 

M Pt = ^M vtx 2 + P* + \P T \. 

where M v f x is the invariant mass of the tracks associated with the reconstructed secondary 
vertex and is the total transverse momentum of the vertex-associated tracks with respect 
to the vertex axis, which we estimated independently of the track momenta by the vector 
along the line joining the IP to the reconstructed vertex position. In this correction, vertexing 
resolution as well as the IP resolution are crucial. 

With these features, topological vertex finding gives excellent 6-tagging efficiency and 
purity. In particular, the efficiency is good even at low P-hadron energies, which is especially 
important for detecting g — * bb. For the selected 4-jet event sample (Section 4) we used our 
simulation to estimate that our mean tagging efficiency for 6-jets in the g — > bb process is 
67%. 



4. Analysis 



Besides the signal events which contain g — » bb, hereafter called 'B events', background events 
can be divided into two categories: 1) events in which a gluon splits to a charm quark pair, 
called 'C events', and 2) events which do not contain any gluon splitting into heavy quarks 
at all, hereafter called 'Q events'. In g — > bb events the two B hadrons from the gluon 
splitting tend to be produced with low energy in a collinear configuration, which allows one 
to discriminate the signal from background. We first required that each event contain 4 jets, 
and that the two jets closest in angle were tagged as b jets on the basis that each contained 
a secondary vertex. We then examined additional kinematic quantities and used a neural 
network technique to improve the signal/background ratio. 



In each event jets were formed by applying the Durham jet-finding algorithm [18| with 



Ucut = 0.005 to the set of charged tracks; this y cu t value was chosen to minimize the sum 
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of the statistical and systematic errors on g h i- Events containing four or more jets were 
retained. The jet-finder was re-run on the > 4— jet events with successively larger y cu t values 
until exactly four jets were reconstructed. With this definition the 4-jet rate in the data was 
(14.58 ± 0.07)%, where the error is statistical only. In the Monte-Carlo simulation the 4-jet 
rate was (14.47 ± 0.02 ± 0.16)% where the first error is statistical and the second is due to 
the uncertainty in the simulation of heavy-quark physics (Section |||). For the B, C and Q 
events the 4-jet rates predicted by the simulation are about 60%, 38% and 14%, respectively. 
Each jet energy was calculated using its associated-track momenta and assuming all tracks 
to have the charged pion mass. 

In each selected event the two jets closest in angle were considered as candidates for 
originating from the gluon splitting process g — > bb, and the topological vertex method was 
applied to them. We required both jets to contain a secondary vertex with a 3D decay 
length greater than 300/im. No tag was applied to the other two jets. 1514 events were 
selected. In each event the tagged jets were labeled '1' and '2', where jet 1 contained the 
vertex with the greater Mp T value, Mp T \, and jet 2 that with the lesser Mp T value, Mp T 2- 
The other two jets in the event were labeled '3' and '4', where jet 3 was more energetic than 
jet 4. With these requirements the selection efficiency for g — > bb events was estimated to be 
16.2%, with a signal/background ratio in the selected sample of approximately 1/10. 75% 
of the background came from Z° — > bb events, 9% from Z° — > qq (g / b) events, and the 
remaining 16% from g — > cc events. In order to improve the signal/background ratio we used 
a neural network technique. We chose the following 9 observables as inputs to the neural 
network; each observable was scaled to correspond to a range between and 1. 

1. Mp T \\ b jets typically have higher values of this quantity than c or uds jets. The 
distribution of Mp T \ is shown in Figure ||. 

2. Mp T 2- This observable has similar discriminating power. The distribution of Mp T 2 is 
shown in Figure ||. 

3. 15Mp T i — P v txi'- where P v t x is the vertex momentum. This observable tends to be 
large for b jets since B decay vertices typically have higher mass than those from charm 
decays, and vertices resulting from B — > D cascade decays have a lower momentum 
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than those from primary D hadrons. The distribution of 15Mp T i — P v txi is shown in 
Figure [| 

4. 15Mp T 2 — P v tx2- This observable also has discriminating power between signal and 
background events. The distribution of 15Mp T 2 — P v tx2 is shown in Figure ||]. 

5. The angle #12 between the vertex axes of jets 1 and 2. The two jets from g — > bb tend 
to have 612 — 0. However, in some Z° — » bb background events a single b jet may be 
split into two by the jet-finder; in these cases the two reconstructed vertices tend to 
have cos #12 > 0.98. The distribution of cos#i2 is shown in Figure ||. 

6. The angle #34 between the axes of jets 3 and 4. In events containing g — > bb this tends 
to be near it, while background events tend to populate the smaller-angle region. The 
distribution of cos#34 is shown in Figure [| 

7. The energy sum of jets 1 and 2, E\ + E2: The two jets arising from g — > bb tend to 
have lower energy than the other two jets in the event. The distribution of E\ + E2 is 
shown in Figure [?]. 

8. The energy sum of jets 3 and 4, £3 + £4: This tends to be larger in signal events than 
in background events. The distribution of £3 + £4 is shown in Figure || 

9. The angle 01234 between the plane 1I12 formed by jets 1 and 2 and the plane H34 



formed by jets 3 and 4: This variable is similar to the Bengtsson-Zerwas angle [19], 
and is useful to separate g — > bb events because the radiated virtual gluon in the 
process Z° — > qqg is polarized in the three-parton plane, and this is reflected in its 
subsequent splitting, by favoring g — > qq emission out of this plane, i.e. 01234 — 7r/2. 



The distribution of COSQ1234 [20] is shown in Figure |9|. 



The measured and simulated distributions agree well for these input observables. We 
trained the neural network using Monte-Carlo samples of about 1800k Z° — > qq events, 1200k 
Z° — > bb events, 780k Z° — > cc events and 50k events containing g — > bb. These samples were 



independent of the ones used for the selection efficiency and background studies. Figure |1C 
shows the distribution of the neural network output variable, Y. We retained events with Y 
greater than 0.7. This value was found to minimise the total error on the final g b i result. 
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Source 


Efficiency (%) 


B 


5.28 ±0.09 (e B ) 


C 


0.165 ±0.018 (e c ) 


Q 


0.00967 ± 0.00038 (e Q ) 



Table 1: Selection efficiencies after all cuts for the three categories. Errors are statis- 
tical only. 

5. Result 

79 events were selected in the data. The number of background events was estimated, us- 
ing the Monte Carlo simulation, to be 41.9, where 35% of the background comes from g 
— > cc events, 63% from Z° — > bb events, and the remaining 2% from Z° — > qq (q ^ b) events. 
Table |] shows the selection efficiencies, relative to the selected hadronic-event sample, for 
the B, C and Q event categories. From these efficiencies and the fraction of events selected 
in the data, /<j = (2.73 ± 0.31) x 10~ 4 , the value of g bb was determined: 

_ fd - (1 - 9cc)eQ - gcc^c ( , \ 

9bb — • K 1 ) 

We obtained 

g bl = (2.44 ±0.59) x 10~ 3 , (2) 
where the error is statistical only. 



6. Systematic Errors 

The efficiencies for the three event categories were evaluated using the Monte-Carlo simula- 
tion. The limitations of the simulation in estimating these efficiencies lead to an uncertainty 
on the result. The error due to limited Monte-Carlo statistics in the efficiency evaluation was 
Ag bi = ±0.12 x 10- 3 . 

A large fraction of events remaining after the selection cuts contain B and D hadrons. 
The uncertainty in the knowledge of the physical processes in the simulation of heavy-flavor 
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±0.24 


b quark mass 


±0.05 


9cc 


±0.12 


Tracking efficiency 


±0.04 


Total 


±0.34 



Table 2: Systematic uncertainties on g b i. 

production and decays constitutes a source of systematic error. All the physical simulation 
parameters were varied within their allowed experimental ranges. In particular, the B and D 
hadron lifetimes, their production rates, and the mean B hadron energy were varied following 



the recommendations of the LEP Heavy Flavour Working Group [21]. We also varied the 
assumed form of the B energy distribution, taking the difference between the optimised 
Bowler and Peterson functions [^] to assign an error. The uncertainties, which are typically 
small, are summarized in Table |2[ 

The dominant systematic error arises via the signal tagging efficiency, eg, from the de- 
pendence on modelling the kinematics of the split gluon: its energy E g , its mass m g and 
the decay angle, 9*, of the two B hadrons in their center-of-mass frame relative to the gluon 
direction. In our default Monte Carlo simulation the kinematics of the signal events are based 
on the JETSET parton shower, which is in good agreement with the theoretical predictions 
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Model g bb (xlO- 3 ) 


JETSET 7.4 


2.20 


PYTHIA 6.136 


2.22 


HERWIG 6.1 


2.50 


ARIADNE 4.08 


2.41 


GRC4F 


2.68 



Table 3: g bb values resulting from different models of the gluon kinematics. 

[]|]. In order to estimate the uncertainty on this assumption, we have investigated alternative 
models, namely the PYTHIA 6.136, HERWIG 6.1 §2j|, and ARIADNE 4.08 @ models, as 



well as the exact leading-order matrix-element calculation GRC4F [25]. In each case, we 
generated, at the parton level, events containing g — » bb. For illustration the E g distributions 



are shown in Fig. 11. The efficiency function computed with JETSET was then reweighted 
by the ratio of the new model to JETSET initial distributions to obtain a new estimate of 
the average efficiency. The resulting g bb values are shown in Table |3| We took the central 
value of this ensemble, g bb = 2.44 xl0~ 3 , as our central result, and assigned a systematic 
error of Ag bb = ±0.24 x 10~ 3 based on the full range of values. For illustration the efficiency 
and its uncertainty are shown in Fig. |ll]. 

The dependence of the efficiency on the 6-quark mass was also investigated at the gen- 
erator level using the GRC4F Monte Carlo program. The variation of e# was computed 
for 6-quark masses between 4.7 and 5.3 GeV/c 2 . The resultant uncertainty is estimated 
to be Ag bb = ±0.05 x 10 -3 . The measured uncertainty in the production fraction of g 
— > cc background events, Ag c5 = ±0.40%, gives an error Ag bb = ±0.12 x 10 -3 . 

In the Monte-Carlo simulation charged tracks used in the topological vertex tag were 
rejected to reproduce better the distributions in the data. Uncertainties in the efficiencies 
due to this rejection were assessed by evaluating the Monte-Carlo efficiencies without the 
rejection algorithm. The difference in the g bb result was taken as a symmetric systematic 
error, Ag bb = ±0.04 x 10~ 3 . 

As cross checks we varied independently the value of y cu t used in the 4-jet event selection 
and the value of the neural network output, Y, used to select the final event sample. In each 
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case we found results consistent with those determined using the optimised value. Table ^ 
summarizes the different sources of systematic error on g bb . The total systematic error was 
estimated to be the sum in quadrature, 0.34 x 1CT 3 . 

7. Summary 

Using the excellent flavor-tagging capabilities of the SLD tracking system, and a new tech- 
nique incorporating a multivariate neural network analysis, we have measured the probability 
for gluon splitting into bb in hadronic Z° decays. The result is 

g bb = (2.44 ± 0.59(stat.) ± 0.34(syst.)) x 1(T 3 

where the first error is statistical and the second is the sum in quadrature of systematic 
effects. This represents the most precise determination of g bb to date. Our result is consistent 
with previous measurements ||. It is also consistent with the theoretical expectation g bb = 
1.75 ± 0.40 x 10 -3 ||], where the central value corresponds to a s {M^) = 0.118 and a 6-quark 
mass of 5.0 GeV/c 2 . Finally, we found that the predictions of the models PYTHIA6.136 (g bb 
= 1.5 x 10~ 3 ) and HERWIG6.1 (g bb = 2.5 x 10~ 3 ), with their default parameter settings, 
are consistent with our measurement. The prediction of ARIADNE4.08 (g bb = 1.3 x 10~ 3 ) is 
slightly below our measurement. 
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Figure 1: Distribution of Mp T \ (see text). The right-most bin includes overflows. The 
points represent the data, the histogram the simulation; the expected contribution 
from g — > bb [g — > cc) (see text) is shown as the dark (shaded) area. 
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Figure 2: Distribution of Mp T 2 (see text). The right-most bin includes overflows. The 
points represent the data, the histogram the simulation; the expected contribution 
from g —> bb [g — > cc) (see text) is shown as the dark (shaded) area. 







+ SLD - 






r~| g-bb E 






J g-cc : 






H-^ □ bkg 


d - 






• 1 











0.2 0.4 0.6 0.8 1.0 

2-2001 

8580A6 (15M Pt1 - P vtx1 )/60 

Figure 3: Distribution of 15Mp Tl — P vtx i (see text). The right-most bin includes 
overflows. The points represent the data, the histogram the simulation; the expected 
contribution from g — > bb [g — > cc) (see text) is shown as the dark (shaded) area. 
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Figure 4: Distribution of 15M Pt2 — P v t X 2 (see text). The right-most bin includes 
overflows. The points represent the data, the histogram the simulation; the expected 
contribution from g — > bb [g — > cc) (see text) is shown as the dark (shaded) area. 
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Figure 5: Distribution of cosOu (see text). The points represent the data, the his- 
togram the simulation; the expected contribution from g — ► bb [g — > cc) (see text) is 
shown as the dark (shaded) area. 
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Figure 6: Distribution of cos#34 (see text). The points represent the data, the histogram 
the simulation; the expected contribution from g —* bb (g — > cc) (see text) is shown as the 
dark (shaded) area. 
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Figure 7: The distribution of the energy sum of jets 1 and 2. The points represent the 
data, the histogram the simulation; the expected contribution from g —>■ bb [g — > cc) 
(see text) is shown as the dark (shaded) area. 
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Figure 8: The distribution of the energy sum of jets 3 and 4. The points represent the 
data, the histogram the simulation; the expected contribution from g —>■ bb [g — > cc) 
(see text) is shown as the dark (shaded) 




Figure 9: The distribution of cos 01234 (see text). The points represent the data, the 
histogram the simulation; the expected contribution from g — > bb (g — > cc) (see text) 
is shown as the dark (shaded) area. 
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Figure 10: Distribution of the neural network output Y. The points represent the 
data, the histogram the simulation; the expected contribution from g — > bb (g — > cc) 
(see text) is shown as the dark (shaded) area. 
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Figure 11: Energy distribution (left-hand scale) of gluons for the g — > bb process in 
different models (see text). Our calculated efficiency and its uncertainty (see text) are 
shown as the band (right-hand scale). 
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